Comment on "Superinsulator and Quantum Synchronization 
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We show that the "theory" of "superinsulating" state presented by V.Vinokur et al (Nature vol. 
452, p. 613, 2008) and Fistul et al (Phys. Rev. Lett. vol. 100, 086805, 2008 ) is essentially incorrect 
due to a sequence of errors in the theoretical analysis of the standard model of Josephson arrays 
which properties have been carefully studied and described in the literature. The line of calculations 
suggested in these articles lead to unphysical results. In particular, the calculations predict a direct 
current flowing through a capacitor. Moreover, this current may flow even in the absence of voltage 
- a sort of supercurrent flowing in the "superinsulating" state. We also question that the theoretical 
model employed in these works, even if treated correctly, is applicable to the analysis of experimental 
data on homogeneously disordered superconductive films. 



1. Two recent articles^* 2 - (hereafter referred to as I 
and II) are supposed to contain theory of "superinsulat- 
ing" state with zero linear resistance first seen in InO^ 
films near superconductor-insulator transition by Sam- 
bandamurthy et al£ and later observed by Baturina et 
al^ in thin films of TiN (unfortunately we found no ref- 
erence in I to the first observation of this phenomenon—). 
Calculations presented in I and II are based on the same 
model and are almost identical to each other both in the 
method and results (few details where they differ will be 
mentioned below), which allows us to analyze them in 
parallel. 

Although the materials studied experimentally in 
Refsi^ are disordered homogeneous films, the authors 
of I and II suggested to model them by a regular 2D 
Josephson junctions (JJ) array. Moreover, they have ar- 
bitrarily chosen a special case when the self-capacitance 
of the islands Co can be neglected with respect to the 
junction capacitance C. This case may be relevant to 
some special geometries of artificial networks but its rel- 
evance to the materials investigated in Rcf. 3 ii has not 
been justified in the papers and is questionable (as will 
be discussed at the end of the present Comment). 

Studying the model of the regular JJ array the authors 
of I, II claim that they discovered a "superinsulating" 
state with a huge dielectric gap (that grows with the size 
of the system) and argue that this state is due to a new 
phenomenon they called "quantum synchronization" . 

In this Comment, we are not going to present our own 
interpretation of the interesting experimental results of 
Refsi^. Our aim here is to demonstrate that the theoret- 
ical part of the paper I, as well as the paper II, contains 
serious errors leading to wrong conclusions. Theoreti- 
cal content of these papers consists both of calculations 
(p. 2 of II and "Methods" part of I) and of a number of 
qualitative arguments. 

We begin with analysis of these calculations and of 
the main result, Eqs. (1) in both papers, then we com- 
ment on some of these qualitative arguments. After that 
we present well-known (for about two decades) results 
on how the model studied in I, II really behaves at low 



temperatures, and finally discuss an applicability of this 
particular JJ array model to the experimental data of 
Refs-iil 

2. The main results of I, II are displayed in the form: 
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and it is stated by the authors of I that Eq. ([I]) goes over 
to Eq. ((2|) at T ^> E c (in fact it does not). Here I is the 
current through the system, V is the voltage and A c is 
the "new" gap estimated by the authors as 

A _ J E c min {A c , L} /d, for ID arrays 
c ~ \ (E c /2) In (min {A c , L} /d) , for 2D arrays > 

and A c = dy/E c0 /E c ^> d is the screening length (E c = 
e 2 /2C, E c0 = e 2 /2Co), L is the sample size; Co is the self- 
capacitance of each island and C 3> Co is the capacitance 
of each junction. 

The "superinsulating" behavior claimed to be found 
in I consists of a very fast (double-exponential) tempera- 
ture decrease of the conductance at low voltage eV < A c 
and low temperatures T < E c , cf. (HJ. Eq. j3]) de- 
scribes a huge (especially for ID case) gap A c depending 
on the size of the system and the authors argue that it 
is formed due to the phenomenon of the "quantum syn- 
chronization" . In fact, as we show below, Eqs. ([U [2|) 
cannot be derived from the suggested model. 

The Lagrangian L of the problem is defined by Eqs. 
(3,4) in I and by Eq. (4) in II. With a small modification, 
keeping only relevant terms considered in I and II, this 
Lagrangian can be written as 



L = L Q + L L 
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and 
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L L = Ej cos (eVt/h + xij (*) + V> (*) ) (6) 
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L « = ^ J2 cos (-eVtyft + XJVj (t) - V (*)) (7) 
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The phases xij and XNj are the leftmost and rightmost 
phases of the sample (next to the leads). 

The authors of I and II write explicitly that, when de- 
riving Eq. ((2|) , they neglect all Josephson couplings inside 
the array, i.e. terms with Ej in L . (cf. paragraphs be- 
tween Eqs. (12) and (13) in I and between Eqs.(8) and 
(9) in II). 

However, once this approximation is made, the system 
becomes electrically disconnected for dc current between 
its left and right terminals. Neglecting the Josephson 
coupling Ej we see that the right and left parts of the 
system are connected to each other by the Coulomb in- 
teraction only (first term in Eq. ([5])), which means that 
the only allowed current is the capacitive one, I(t) cx 
C c gdV(t) I 'dt, which is strictly zero in the dc limit. On a 
more formal level, the neglect of all the intrinsic Joseph- 
son couplings leads to the possibility of rotating the phase 
difference xij W — XNj(t) by an arbitrary constant with- 
out any energy cost (gauge invariance). Writing the 
corresponding expression for the current and integrating 
over static \in or XNj in the interval (0, 2ir) one comes 
to the strictly zero dc current. 

How could it happen that the authors of I, II got the 
nonzero direct current? We believe that it is due to an 
incorrect approximation made during the derivation of 
Eqs. (11) and (7) of the articles I and II, respectively. 
Namely, they first represent the sum of two cosine terms 
Lr + Ll in the equivalent form of the product (and then 
the same with the expression for the current): 



M 

2Ej ^2 cos 



Xi ,j + XN,j 



(9) 



2eVt + ip + xij - XNj 



and then they replace the average of the product of two 
nontrivial functions of time (e.g. Eq.(7) of I) by the prod- 
uct of their averages, thus coming to Eq. (11) of I (or 
Eq.(7) of II). We rewrite these equations for the current 
I s (V) here for an easy reference 



I s (y) = Aim / dt exp [St + 2ieVt] K (t) 
Jo 



(10) 



where A is a constant proportional to Ej, and the func- 
tion K (t) equals 

K (t) = (exp i [xi (t) - xi (0) - X n (t) + Xn (0)]), 

(11) 

and t is real time. 

This procedure of I and II could make sense under some 
special circumstances. Indeed, in the situation studied in 
Refs. [18,19] of I, there are no intrinsic array junctions 
and the autocorrelation function K(t) = 1. The same 
approximation could also be used if the Josephson cou- 
plings of the bulk were different from those near the leads 
and were sufficiently large (i.e. if the condition Ej ^> E c 
would be fulfilled in the bulk of the array). In this case, 
phase fluctuations in the bulk would be strongly sup- 
pressed and this situation would not be essentially differ- 
ent from the one considered in Refs. [18,19]. This means 
that Eq. (TTU)) can in some cases (when K (t) is slow func- 
tion of time) be correct and the calculation based on the 
analog of the representation (1. 11) is, in principle, possi- 
ble. 

However, an attempt of Vinokur et al (I) and Fistul 
et al (II) of using formulae written in [18,19] to their 
own model failed explicitly due to the neglect of cross- 
correlations between two trigonometric factors contained 
in Eq.7 of I for vanishing values of the Josephson cou- 
plings in the bulk. The function K (t) calculated in I, 
II changes with time extremely fast (on the scale of 
A^ 1 ) and there is no justification of the "decoupling" 
approximation employed. The correct calculation within 
their approximation of zero Josephson couplings inside 
the sample would lead the authors to the evident and 
trivial result: dc current cannot flow through a capaci- 
tor. 

Furthermore, Eqs. (|Tl [3J) were written in I, II without 
specifying the form of the pre-exponential. Following the 
calculation of I, II one comes to the conclusion that the 
pre-exponential is proportional to Ej (see Eqs. ((TOlfTT]) ). 
This is incorrect: expanding in the coupling terms of Lq 
one can see that the coefficient in front of such a term 
is strictly zero if N > 4. In fact, using the initial La- 
grangian, Eqs. (0HH]), a non-zero contribution to the cur- 
rent could be obtained (neglecting all inelastic processes 
as the authors of I and II did) in the JV-th order of the 
expansion in a small parameter Ej /E c only (where N is 
the size of the array). 

One could wonder, however, if Eqs. (1,2) could still 
be derived, although with a different coefficient (not just 
proportional to Ej). 

We found that Eqs. fl] ^ cannot be repaired in such 
a way. The form of Eqs. JTJ) and ((2|) leads to the expo- 
nential suppression of the current not only at low voltage 
eV < A c but at eV > A c as well. This surprising phe- 
nomenon was not commented in paper II. At the same 
time, in the later paper I the authors mentioned that 
they used their result (1) for eV < A c only. However, 
such a limitation does not follow from any step made in 
their derivation of this result and the need to invoke it 
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on a purely verbal level indicates erroneous nature of the 
analysis. 

Another problem occurs when the authors neglect con- 
tributions of non-zero winding numbers at T larger than 
E c in Eq. (1.13). The quantization of charge is inherently 
related to the summation over all winding numbers. No 
insulating behavior can be obtained without it. 

Actually, the necessity of taking into account the dis- 
creteness of the charges even for T > E c is specific for 
the model with the vanishing self-capacitance Co consid- 
ered in the articles I, II. Adding the elementary charge 
2e to the system costs the huge energy A c , Eq. ([3]), 
and not E c . Therefore, one must sum over the integer 
charges for all temperatures T < A c and the substitu- 
tion of the sum by an integral is impossible. In terms of 
actual calculations it means that the integral over non- 
zero Fourier-components \u> must be complemented by 
the sum over nonzero integer winding numbers. 

Neglecting winding numbers the authors take the 
Gaussian integral over the phase \ an d obtain (see Eq. 
(10) of II). 

K (t) = exp (-2A C E C ( (T) t 2 - 2iA c t) (12) 

where ( (T) = T/E c for T > E c . 

Note the presence of t 2 term in Eq. (|12p. It is this 
term that might lead to Gaussian integrals and, eventu- 
ally, to Eqs. ([I] [2]) for some physical quantities. However, 
the summation over the non-zero winding numbers k re- 
sults in the vanishing of the t 2 term in the exponent in 
Eq. (fT2|) . This is explicitly demonstrated, e.g., in Ref. 5 , 
where a granular metal consisting of normal grains was 
considered. In the absence of the t 2 term in the expo- 
nent in the function K (t) one would not be able to get 
generally the combination (A c — eV) 2 in the exponents 
in Eqs. (UJUD. 

There is an exception when the correlation function 
K (t) can really be written in a Gaussian form for T < 
A c . This is the case of a one dimensional chain of J J with 
the vanishing self-capacitance Co. In this case, the corre- 
lation function K(t), Eq. (fTTj) . factorizes into a prod- 
uct of correlation functions for each junction between 
the grains. The authors of I, II use this property and 
bring the expression for the correlation function K (t) 
to the form of Eq. (|12p for arbitrary temperature with 
a function £ (T) proportional to T at T > E c , and to 
exp (— E c /T) at T < E c . Then, apparently calculating 
Gaussian integrals they come to Eqs. |TJ[2]), thus obtain- 
ing the double-exponential behavior at low temperatures. 
The latter is supposed to describe the "superinsulator" . 

However, direct inspection shows that Eqs. tfTOl [T2| do 
not lead to the results announced in Eqs. ([HE]) (&t this 
moment we do not even discuss the validity of (fTPI [T2")) ). 

Substituting Eqs. {TJD into Eq. (JTOJ) (Eqs. (14) or (15) 
of I into Eq. (11) of I) we obtain 

I a (V)=A / sin [2 (V- A c ) t] exp (St) K (t) dt, (13) 
Jo 



where Kq (t) is a real, even function of time defined via 
the relation K (t) = exp (-iNE c t/2) K (t). 
We see immediately from Eq. (fTBl . that 

I. (V) ± -I. (-V) (14) 

In particular, Eq.fTJJ shows that the "superinsulating" 
state can carry a non-zero current without any voltage 
or magnetic field applied, i.e. the ground state maintains 
a "supercurrent" : 

/>OC 

I s (0) = -A sin [(2A c t)] exp (St) K (t) dt (15) 
Jo 

It is not difficult to see that this supercurrent is even not 
small. Using Eq. (fl"2|) we obtain in the limit £ ^ 

I. (0) ~ -A/ A c (16) 

The only way to obtain Eqs. HI [3J) from Eqs. (Q31 [12]) 
is to use ImK (t) instead of K (t) in Eq. (fTU|) . But this 
is hard to justify. After such a replacement, for the two 
junction system, where K (t) = 1 (in this case Eq. (TlT)]) 
is valid and gives a non-zero current consistent with pre- 
vious works), one would get just zero current after the 
substitution K — > lmK(t). The authors of I, II em- 
phasize explicitly (between Eqs. (12, 13) in I and Eqs. 
(8, 9) in II) that Eq. ((TOj) written with K (t) (not with 
ImK (t)) allows them to recover the results for the two 
junction systems. 

So, we see that following the calculational scheme of 
I, II one comes inevitably to the unphysical results, Eqs. 
(IT31 [TB|) . In fact, this is a result of one more mistake 
made by the authors. It is clear that the non-zero "su- 
percurrents" in the "superinsulator" originate from the 
non-zero value of ImK (t) obtained for the real time t. 
However, it follows from the definition of the function 
K(t), Eq. (fTTj) . that its imaginary part is identically 
zero because the Lagrangian Lq is an even function of x, 
whereas the imaginary part of the exponential entering 
Eq. (|11|) is odd in \- Any type of averaging with Lq 
gives in this situation zero for ImK (t). 

Actually, the non-zero value of Im K (t) obtained in 
I, II is a consequence of an erroneous analytical con- 
tinuation from the imaginary time r to the real one t. 
The standard procedure of the analytical continuation 
for computation of time-dependent quantities is as fol- 
lows: i) to calculate the proper correlation function at 
Matsubara frequencies, ii) to make analytical continua- 
tion to real frequencies, and then, iii) to perform Fourier 
transformation to real times. Simple replacing r — > it, as 
was done in I, II, is not correct procedure. 

Now we review few of "qualitative arguments" given 
in I and II. 

1) In the last paragraph of left column of p. 3 of II, as well 
as in various locations in I, the authors write that they 
obtained "global phase-synchronized state..." after they 
used the model containing no Josephson couplings inside 
the array. On the other hand, the authors are aware of 
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the duality relation between the charge and phase vari- 
ables (cf. left column of p. 3 in paper I), which means the 
following statement: the existence of strong phase corre- 
lations implies weakness of charge number correlations, 
and vice verse. We are not able to understand the log- 
ics of applying (simultaneously) both the concepts of the 
charge quantization (which is the essence of the Coulomb 
blockade) and the "synchronization of phases" . 

2) On the first page of I we see an estimate E c ~ A/g, 
where g > 1 is a dimensionless intergrain conductance, 
i.e. E c is smaller than the superconducting gap A and 
thus smaller than the Josephson coupling Ej ~ <?A. We 
found no obvious arguments that would allow one to ob- 
tain this estimate in a model which explicitly assumes 
E c > Ej, as we see in "Methods" part of I. 

3) The statement in the 1st paragraph of the right col- 
umn, p. 3 of II reads: "... even large... fluctuations in 
E c , E c o and Ej as well as offset charges are negligible 
as compared to the huge magnitude of A c " . We explain 
below that if any threshold voltage Vr — A c /e exists, it 
must scale with the length of the system, irrespective of 
any model details. Thus, the statement cited above con- 
tains a comparison of intensive (size-independent) quan- 
tities with extensive ones. To give an example: if such 
a statement could be correct, any local disorder would 
be irrelevant to the thermodynamics properties of any 
macroscopic system. 

4) The first paragraph of p. 4 in paper II contains a com- 
parison between the lowest-temperature threshold volt- 
age Vt and activation gap To measured ins at higher tem- 
peratures. In order to "explain" the high ratio eVr/T 
the authors invoke an idea of "dielectric breakdown" and 
use for 2D system results of calculations obtained for ID 
model. They did it since for their 2D model they obtained 
threshold voltage Vr — A c /e that scales just logarithmi- 
cally with L, cf. second line of ([3]), which is not enough 
to get the very large ratio eVr/T . 

We note that: i) nothing in the model studied by 
the authors indicates the phenomenon of the "dielectric 
breakdown" raised by the authors, ii) there is no reason 
to invoke such a "concept" since the scaling of the thresh- 
old voltage with the length L is a direct consequence of 
any correct calculation leading to a nonzero Vr- The rea- 
son for that is simple: the relevant intrinsic parameter for 
a macroscopic system is the electric field £ and a criti- 
cal value £t of this field determines the threshold. The 
value of £t does not depend on L for large L, whereas 
the voltage Vr oc L. The authors of I, II failed to obtain 
such scaling in 2D due to several calculational mistakes 
analyzed above. 

5) At last, we mention an explanation given in II, cf. Eq. 
(14), for the non-monotonic dependence of the insulat- 
ing gap on magnetic field: it is based on a replacement of 
the disordered superconducting media by a single (!) dc 
SQUID, with its periodic dependence on magnetic flux 
through the SQUID loop. Remarkably, this "explana- 
tion" disappears completely from paper I, being replaced 
there by the statement that the magnetic field suppresses 



the "superinsulating state" like it does with the supercon- 
ducting state. That version sounds better in view of the 
experimental data presented in Fig. 2, but it has no rela- 
tion to the calculations based on the JJ array model. In 
addition, it does not give any hint on the non-monotonic 
independence observed ini^ and "explained" in II. 

3. Now we remind the readers some known properties 
of the model of I, II derived in a number of earlier pa- 
per^™— — , cf. also a review article^: 
a) The insulating regime is realized at low temperatures 
if Ej < aE c , where a ~ 1. In this case the adequate 
description of charge transport involves rare hopping of 
Cooper pairs between the islands, as was explained al- 
ready in Ref£. 

The number of these charge excitations is exponen- 
tially low at T <C E c , leading to the exponential sup- 
pression of the conductivity. The Coulomb interaction 
between the Cooper pairs in the grains was described 
in Ref. 6 using a general capacitance matrix dj without 
making any assumption about its parameters. This is 
why it was assumed that the energy of adding a particle 
into the system was of the order E c . 

b) The model with the vanishing self-capacitance, 
Co = 0, employed by the authors of I, II is special because 
the energy of interacting charges 

where p% are charge operators for the Cooper pairs, is at 
large scales linear in ID and logarithmic in 2D. This is 
precisely at the origin of the "big gap" A c found in I, 
II and displayed in the original form in Eq. Q. Note, 
however, a wrong coefficient in the second line for 2D: 
in fact, the energy of a single 2e excitation is equal 
to ~E c ln(L/d). The logarithmic interaction of charges 
leads to a charge-binding transition of the Berezinsky- 
Kosterlitz-Thouless (BKT) type. The temperature T c of 
this transition was estimated in Refsi^ ancU^. Upon 
approaching this temperature from above, the linear re- 
sistivity diverges, cf. e.g. Eq. (9) of RefJ£. 
c) At T < T c all the charge excitations are paired and 
linear conductivity vanishes, in "dual analogy" § with the 
vanishing of the linear resistivity in 2D superconductors 
below BKT transition. Nonlinear I(V) transport at low 
voltage is due to electric-field induced unbinding of neu- 
tral "charge molecules" , 

I(V) oc (£/£ T ) a a oc T c /T £ = V/L, (18) 

again similar to V(I) characteristics in 2D superconduc- 
tors below BKT transition; the threshold electric field 
can be estimated as £t ~ E c /ed. The corresponding 
characteristic voltage Vt = L£t scales with L for triv- 
ial reason, being extensive variable (just like the total 
current through a superconductive 2D array scales with 
its width L±), without any "quantum synchronization" 
"discovered" in paper I. 
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4. We have arrived at the conclusion that the results 
obtained in I, II and expressed by Eqs. flj [2]) are in- 
correct. One may wonder if, however, some qualitative 
features of them might still be correct. Of a particular 
interest is the "superinsulation" formula represented by 
the low-V limit of Eq.©: 



R tx exp 



A, 



cxp(£ c /2T) 



(19) 



Since this formula was suggested in paper I for ID case 
only, we analyze this specific case below and show that 
Eq. (|19p is incompatible with simple physical arguments. 

In ID case A c oc L, and thus, Eq. p^]) demonstrates 
both the exponential dependence on the system size and 
much faster than exponential dependence on tempera- 
ture. If such a dependence on T and L could indeed be 
obtained, it would be really unexpected. We show now 
that such a behavior does not exist. 

The ground state energy is achieved when there are no 
charges in the system. Then Eq. (fl"7|) gives zero Coulomb 
energy U — 0. A configuration with the lowest non-zero 
energy is a dipole consisting of the charges +2e and — 2e 
located at neighboring grains. In the model considered 
here the energy U (x) of the dipole of the length x linearly 
depends on x and, subjected to the external electric field 
£, equals 



U{x) = 2e(£ c - \£\)x 



(20) 



The energy of the dipole linearly grows with x for \£\ < 
£ c . A current through the system is possible provided 
the size of the dipole can reach the size of the system L. 
The energy corresponding to such a dipole is huge: 



U (L) = 2 (A c - eV) = 2eL(£ c - \£\) 



(21) 



where £ c — A c /eL. However, the probability of creat- 
ing this dipole is finite at finite L and proportional to 
cxp(—U(L)/T), therefore the resistivity R can be writ- 
ten as 



R oc exp 



2eL (£ c -\£\) 
T 



(22) 



We see from Eq. (|22j) that the resistivity obtained 
from this simple consideration contains the combination 
A c — eV and not square of it. Therefore, the Coulomb 
blockade is important for eV < A c only. 

Eq. (|22p shows that the Coulomb blockade can be 
overcome if one creates a dipole with the energy of order 
A c . Clearly, this probability is small but it is described 
by the activation law. This simple argument excludes 
the double-exponential temperature dependence of the 
resistivity R, Eq. 

We want to emphasize that the big gap in the spec- 
trum of the excitations A c exists (within the specific ID 
model discussed) not only for the Cooper pairs but for 
the normal electrons as well^ (in the latter case one would 



have to use A c /4 instead of A c ). So, if the "superinsu- 
lation" was due to an arbitrarily chosen special form of 
the capacitance matrix, it would exist even without any 
superconducting pairing. 

An important remark is in order: following I and II we 
considered the situation of an exponential dependence 
of resistivity on the system size L. In reality, this usu- 
ally does not make sense for large systems: in 2D case 
a resistance per square p{T) can be defined, whereas in 
ID case the resistance grows oc L. This behavior is due 
to inelastic processes including thermal activation. Such 
processes (neglected within the calculational scheme used 
in papers I and II) lead, in particular, to the current- 
voltage dependence (fT5|) for 2D case. ID model with 
weak Ej and vanishing self-capacitance is very specific 
due to its simplicity: it reduces to a single capacitor with 
C c fF = C/N. Thus, the meaning of Eq. (|2"2"j) is also triv- 
ial: it corresponds to a resistance through a capacitor 
in the Coulomb blockade regime. Actually, non-trivial 
macroscopic systems with a low-temperature conductiv- 
ity proportional to exp (— cL) do exist: such a behavior 
is a consequence of a "quantum topological order" , cf. 
e.gJ^. The standard 2D Josephson-junction model dis- 
cussed in I, II does not lead to such a behavior. 

5. Could this 2D regular J J array model be used to de- 
scribe experimental data&i on homogeneously disordered 
superconductive films? The authors of I wrote (top of 
the last column of the main text) : "our understanding of 
the origin of superinsulating state in the films relies on 
the formation of the network of superconducting droplets 
within the normal matrix" . Still, the model they consid- 
ered contains superconductive islands connected by tun- 
nel barriers without any sign of normal matrix inside. 
The problem of the superconductive islands put into a 
normal matrix is essentially different in its physics and 



13.14,15,16 



m corn- 



more complicated for theoretical studies 
parison to the model of 2D JJ array; its relevance to the 
data of^ could be a subject of separate work. Coming 
back to the discussed model of 2D JJ array, we note its 
serious deficiencies regarding a possibility of an applica- 
tion to homogeneously disordered superconductors: 

i) in the absence of well-defined (structural) grains sep- 
arated by tunnel junctions, the notion of the charging 
energy E c becomes ill-defined, to begin with; 

ii) even if the grains separated by tunnel barriers are as- 
sumed to appear "somehow" in a continuously disordered 
media, there is no reason to expect them to form a reg- 
ular lattice with identical areas, Ej and E c . 

iii) even if in spite i) and ii) a granular model could with 
some degree of imagination be thought of as relevant, the 
model considered in I, II contains a vanishing ratio of the 
capacitances Cq/C and this is the crucial point to obtain 
the "huge" gap; why should this condition be fulfilled in 
a "self-organized" granular network? 

iv) care should be exercised with respect to the role of 
random stray charges, which are always present in real 
arrays and influence strongly the insulating part of the 
phase diagram for any model where the insulating be- 
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havior is due to Coulomb repulsion (cf. e.gi 17 ' 18 and 
references therein). In particular, in the presence of a 
charge disorder with a typical amplitude of order e, the 
charge-pairing BKT-like transition does not exist. 

Based on these arguments we question any kind of "fit- 
ting" the experimental data — by results of any calcula- 
tions based on the simplified JJ array model. 



To summarize: the phenomenon dubbed as "quan- 
tum synchronization" does not exist in the model stud- 
ied. This model describes the well known phenomena of 
Coulomb blockade. The result obtained by the authors 
are in conflict with previous studies of the model and we 
believe are based on erroneous calculations. 



1 V. M. Vinokur, T. I. Baturina, M. V. Fistul, A. Yu. 
Mirononv, M.R. Baklanov,and C. Strunk, Nature, 452, 613 
(2008) 

2 M. V. Fistul, V. M. Vinokur and T. I. Baturina, Phys. 
Rev. Lett. 100, 086805 (2008). 

3 G. Sambandamurthy, L.W. Engel, A. Johansson, E. Peled 
and D. Shahar, Phys. Rev. Lett. 94, 017003 (2005). 

4 T. I. Baturina, A. Mironov, V. M. Vinokur, M.R. Bak- 
lanov, and C. Strunk, Phys. Rev. Lett. 99, 257003 (2007). 

5 K.B. Efetov and A. Tschersich, Phys. Rev. B 67, 174205 
(2003) 

6 K.B. Efetov, Zh. Eksp. Teor. Fiz. 78, 2017 (1980) (Sov. 
Phys. JETP 51, 1015 (1980)) 

7 J. E. Mooij et al, Phys. Rev. Lett. 65, 645 (1990) 

8 M. P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 

9 R. Fazio and G. Schoen, Phys. Rev. B 43 5307 (1991) 

10 M. V. Feigel'man, S. E. Korshunov and A. B. Pugachev 
Pisma v ZhETF, 65, 541 (1997) [JETP Lett. 65, 566 
(1997)]. 

11 R. Fazio and H. van der Zant, Phys. Rep. 355, 235 (2001). 

12 L.B.Ioffe and M.V. Feigel'man, Phys. Rev. B 66, 224503 
(2002). 

13 M. V. Feigelman and A. I. Larkin, Chem. Phys. 235, 107 
(1998) 

14 B. Spivak, A. Zyuzin and M. Hruska, Phys. Rev. B 64, 
132502 (2001) 

15 M. V. Feigelman, A. I. Larkin and M. A. Skvortsov, Phys. 
Rev. Lett. 86, 1869 (2001). 

16 M. A. Skvortsov and M. V. Feigelman, Phys. Rev. Lett. 
95, 057002 (2005) 

17 E. Granato and J. M. Kosterlitz, Phys. Rev. B 33, 6533 
(1986). 

18 I. V. Protopopov and M. V. Feigelman, Pisma v ZhETF 
85, 621 (2007) [JETP Lett. 85, 513 (2007)]. 



I. NOTE ADDED ON JULY 30, 2008. 

Recently, M. V. Fistul, V. M. Vinokur, T. I. Batu- 
rina (FVB) submitted a more detailed paper III (Macro- 
scopic Coulomb blockade in large Josephson junction ar- 
rays, arXiv:0806.4311, vl). Several formulae criticized 
in our comment have been corrected. In III the current 
I s (V) is expressed through a retarded correlation func- 
tion of hopping amplitudes and it is shown how to write 
this function at real time. This makes the current an odd 
function of the voltage V as it should be. 

However, the most serious errors remain. A major ap- 
proximation when the first cosine in Eq. (9) (of our Com- 



ment) is replaced by a constant would be justified if the 
Josephson couplings in the bulk were large. However, 
the authors consider an opposite limit. There the phase 
difference between grains in the bulk strongly fluctuates 
with time and the stationary perturbation theory used in 
III, Eqs. (10-13), is not applicable. For the same reason, 
the arguments lead to Eq. (31) of III are not consistent. 
Setting all couplings in the bulk to zero, as the authors 
do, breaks contacts between grains. This system (thus 
approximated by a capacitor) cannot support a finite dc 
current. 

Nevertheless, in papers I-III the authors present for- 
mulae for dc current as a smooth function of voltage V 
(see Eqs. (46, 54) of III). How could it happen? 

We have found that an integral obtained as a result of 
their approximation has not been calculated properly. 

Following I-III one can write the current in the form 
of Eq. (25) of III with the correlation function K (t) 
written in the first equation of Eqs. (33) of III (or, Eq. 
(52) of III). For a large number of grains one can neglect 
an influence of the leads and put Ki ea d s (t) = 1. The 
integral obtained in this way can be easily evaluated in 
different limits. 

For example, in the limit E c /lnN <C T <C E c , where 
N is the number of the grains in the system, the authors 
of III obtained their Eq. (53). They calculate the inte- 
gral using the saddle point method: expanding cos (E c t) 
entering Eq. (53) in t near t = and then computing 
the Gaussian integral around this saddle point. In this 
way the authors obtained the double exponential behav- 
ior (see Eq. (1) of our comment). 

However, t = is not the only saddle point of this in- 
tegral: all t — ±27rm with integer m are saddle points as 
well, and all of them must be taken into account. Inte- 
grating near all the saddle points and summing them up 
yields the following expression for the current 



Is(V) 



E C /4T 



N 



(Jo (V) - J (-V)) 



(23) 



where 



00 / k 2 \ 

h[V)= S(2eV~E c (N/2 + k))e Xp (-—e E ^ T ) 



and B is a coefficient proportional to Ej. 



Eq. (|23|) contains a sum of matrix elements between 
the states with different numbers of the particles of iso- 
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lated grains. The discreteness of the levels is the con- 
sequence of the discreteness of the charge of the Cooper 
pairs. It is there at all temperature regimes once contacts 
are broken. The levels may be smeared only by taking 
into account charge tunnelling between the grains which 
has been explicitly ignored by the approximation made 
by the authors. A tunnelling (a non-zero Ej) smoothes 
/(V) functions, not the temperature. This, seems to be 
an obvious point, has been discussed in the main part 
of our Comment (see e.g., Eqs. (18,22)). No disorder 
can help in this respect as long as the contacts remain 
broken. 

Correct calculation of the integral does not save Eq. 
(|23|) . though. It still shows a dc-current flowing through 



a capacitor at resonant values of voltage. 

Finally we comment on the last paragraph on p. 11 of 
III. There the authors explain how one should think of 
"superinsulating" state: the internal part of the array 
(except the rightmost and leftmost junctions) "acts co- 
herently as a single superconducting island" . It is not 
quite clear whether the authors have in mind a closed 
ring to be a superconductor or "superinsulator" . Re- 
gardless of this, we do not think that superconductivity 
in an array with broken junctions is possible anyway. 

We conclude that the new theoretical material pre- 
sented in III neither clarifies nor essentially corrects the 
analysis and the results of I and II. 



